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Nonomura and Hu Reply: In the preceding comment 
, Olsson and Teitel questioned the possible vortex slush 
(VS) phase in the frustrated XY model with point defects 
reported by the present authors 0- The VS phase was 
originally proposed in order to explain an experiment of 
irradiated YBazCuaOy-,? (YBCO) 0. This phase was 
also observed in an optimally-doped pristine YBCO 0, 
where the VS phase locates above the Bragg glass (BrG) 
phase in the H-T phase diagram. In Monte Carlo simula- 
tions of the frustrated XY model by the present authors 
, a first-order transition was observed between the vor- 
tex liquid (VL) and VS phases up to a certain density of 
point defects. The structure factor in the VS phase shows 
obscure Bragg peaks, which was interpreted as a short- 
range order in the ab plane. In comparison with the BrG 
phase for lower density of point defects, the VS phase 
shows much larger density of dislocations in the ab plane 
and the vanishing helicity modulus along the z axis. 

Olsson and Teitel simulated the same model for the 
same parameterization as used in Ref. Q. On the ba- 
sis of the structure factors observed in layer by layer, 
they argued that the VS phase observed by the present 
authors may be an artifact of a finite system size, and 
that this region may be included in the BrG phase in 
the thermodynamic limit. However, their argument in 
the Comment is not justified sufficiently by the provided 
numerical data for the following reasons. 

First, they observed strong hysteresis behavior by 
sweeping the pinning strength of point defects e in the 
VS region, and took this behavior as the evidence for a 
wide coexistence region of the first-order melting of the 
Bragg glass. However, this behavior can alternatively 
be interpreted as merging of the consequent VL-VS and 
VS-BrG first-order phase transitions due to a small sys- 
tem size. The two-step behavior in the hysteresis curve 
of the peak value of the structure factor in Fig. 1 (b) of 
Ref. looks consistent with the latter picture. It should 
also be pointed out that the hysteresis behavior may be 
enhanced by their e-sweeping procedure. Since the an- 
nealing process is not included in this procedure in spite 
of possible drastic changes of the configurations of flux 
lines caused by varying e, Monte Carlo steps necessary 
for equilibration in this procedure may be much larger 
than those in the temperature sweeping adopted in our 
previous article |2J. 

Second, they argued that the energy loss due to a mis- 
match in different layers is proportional to J Z L 2 with the 
transverse system size L. This scaling argument is based 
on the assumption that the mismatch characterized by 
the change of peak positions of Bragg peaks occurs as 
abruptly as a domain wall of the Ising model. How- 
ever, the mismatches displayed in Fig. 2 of Ref. Q relax 
across a number of layers. When the relaxation takes 
place across L w layers, the energy loss is proportional to 
J Z L 2 /L W . It is natural to expect that L w depends on 
the thickness of the system L z . Provided L w is propor- 


tional to L z ~ L, the energy loss caused by a mismatch 
is proportional to J Z L, instead of J Z L 2 . Then, their con- 
clusion should be changed completely ||. In order to 
address the issue sufficiently, one should vary the sys- 
tem size and check the size dependence of the relaxation 
between mismatches and the number of mismatches. 

On the other hand, we have to say that our previous 
study cannot completely exclude the possibility that 
the "first-order VL-VS phase boundary" may actually 
be the melting line of the finger-like wiggled BrG region 
stretching into the VL phase 0. Such a narrow BrG 
region (see Fig. 1 of Ref. 0) may be weakened by ther- 
mal fluctuations and the inverse-melting behavior may be 
masked by finite-size effects in numerical simulations or 
by some experimental conditions. Even if it is the case, 
the shape of the vortex phase diagram is qualitatively 
different from that obtained by Olsson and Teitel pj. 

We also notice that the vortex phase diagram including 
the VS phase as Ref. p| was also reported in a YBCO thin 
film 181 , and that there exist several theoretical studies 
pl llfilll[ rL2j| consistent with our numerical results. 

In summary, the scaling argument which plays a key 
role in the preceding Comment cannot be sufficiently jus- 
tified by the provided numerical data. Their data are 
indeed consistent with our previous article. 
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